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Beta decays of the isotones with N=126 are studied by shell model calculations taking into account 
both the Gamow- Teller (GT) and first-forbidden (FF) transitions. The FF transitions are found 
to be important to reduce the half-lives, by nearly twice to several times, from those by the GT 
contributions only. Possible implications of the short half-lives of the waiting point nuclei on the 
r-process nucleosynthesis during the supernova explosions are discussed. A slight shift of the third 
peak of the element abundances in the r-process toward higher mass region is found. 
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I. INTRODUCTION 



The r-process is the most promising process for the 
synthesis of about a half of heavy elements beyond 
iron Q. Study of the r-process element synthesis has 
been done by considering neutrino-driven winds in su- 
pernova explosions as well as ONeMg supernovae Q 
and neutron-star meregers The production of the 
r-process elements in collapsars relating to long gamma- 
ray bursts has also been investigated (sj . Recent observa- 
tions of r-process elements in extremely metal-poor stars 
have suggested that the r-process has occured in ex- 
plosions of massive stars or in neutron-star mergers from 
early stage of Galactic chemical evolution. The r-process 
is affected by various inputs of nuclear properties such 
as mass formulae Q, /3-decay rates, (n, 7) and (7, n) 
reaction cross sections, reaction rates of a-processes in 
light nuclei [1], etc., and depends also on astrophysical 
conditions; electron-to-baryon number ratio, Y^, the en- 
tropy and temperature of the explosion environment, and 
neutrino processes d, [l3|- The evaluation of /3-decay 
rates, particularly at the waiting point nuclei, is one of 
the important issues of the nucleosynthesis through the 
r-process. Investigations on the ^-decays of isotones with 
neutron magic number of N=82 have been done by var- 
ious methods including shell model [11], QRPA/FRDM 
[11], QRPA/ETFSI lii and HFB-hQRPA ^ calcula- 
tions as well as CQRPA ones [l^. The half-lives of nu- 
clei obtained by these calculations are rather consistent 
to one another, and especially in shell model calculations 
experimental half-lives at proton number of Z =47, 48 
and 49 are well reproduced [ll|. 

For the /3-decays at A^=126 isotones, however, half- 
lives obtained by various calculations differ to one an- 



other [13] ■ First-forbidden (FF) transitions become im- 
portant for these nuclei in addition to the Gamow- Teller 
(GT) transitions in contrast to the case of N—82. A 
strong suppression of the half-lives for the A^=126 iso- 
tones due to the FF transitions has been predicted in Ref. 
[ist . Shell model calculations of the /3-decays of 7V=126 
isotones have been done only with the contributions from 
the GT transitions [13, [3 • Moreover, experimental data 
for the /3-decays in this region of nuclei are lacking. The 
region near the waiting point nuclei at 7V=126 is, there- 
fore, called 'blank spot' region. 



Here, we study /3-decays of iV=126 isotones by taking 
into account both the GT and FF transitions to evaluate 
their half-lives. Shell model calculations are done with 
the use of shell model interactions, modified G-matrix el- 
ements, that reproduce well the observed energy levels of 
the isotones with a few (2 to 5) proton holes outside ^ospj^ 
[13, [II- In sect. II, GT and FF transition strengths 
for the A^=126 isotones are evaluated, and half-lives of 
the /3-decays are obtained by taking into account both 
the transitions. The calculated half-lives of the waiting 
point nuclei in the r-process become shorter by including 
the contributions from the FF transitions. In Sect. Ill, 
possible implications of these half-lives for the r-process 
nucleosynthesis are discussed. A summary is given in 
Sect. IV. 
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II. TRANSITION STRENGTHS AND 
HALF-LIVES OF THE iV =126 ISOTONES 

A. Gamow- Teller, Spin-dipole Transition Strengths 
and Shape Factors 

First, the method for the evaluation of the rates and 
half-hves of the GT and FF transitions are explained. 
Then isotone dependence of the GT and FF transition 
strengths as well as the effects of the inclusion of the FF 
transitions on the half-lives of the isotones are discussed. 

The decay rates, A, as well as the partial half-lives, 
ti/2, of the transitions are obtained by the following for- 
mulae fiiSI; 



A = ln2/ti/2 = //8896 (s^i) 



C {w)F {Z , w)pw{'Wo — vSf'dw 



C{w) = Ka + Kiw + K-i/w + K2W^ 



(1) 



where w is the electron energy, F{Z,w) is the Fermi func- 
tion, and KnS depend on nuclear transition matrix ele- 
ments. Here, relativistic corrections from the expansion 
of electron radial wave functions in powers of electron 
mass and nuclear charge parameters are included; ma- 
trix elements of one-body operators, [r x t?]^ with A =0, 
1, 2 and r, as well as those from weak hadronic currents, 
75, a, are taken into account for the FF transitions. In 
case of the GT transition, the shape factor C{w) does 
not depend on electron energy; 



1 



2J, + 1 
0(1+) = gACTt- 



</l|0(l+)|N> 



(2) 



with gA the axial vector coupling constant, Ji is the spin 
of the initial state and t^\n > =\p >, and {n ^0) —0. 
In case of the FF transitions, for 0~ transitions. 



distribution on the electron wave function are taken into 
account 111, u3. For 1~ transitions. 
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(4) 



where is the nucleon mass and the primes in x' and 
u' indicate that the effects of finite nuclear charge distri- 
bution are taken into account [l^ [l^ . The quantities 
and A2 are defined "2^ in terms of electron wave functions 
and depend on electron momentum, p^. Their values are 
close to unity; ^1 =0.9^1.0 and A2 = 0.7^1.0 at Pe > 
0.5 and gets larger than 1 at pe < 0.5 for Z considered 
here 

For 2~ transitions. 



K2 = -z\l + X2) 

z^2gA< f\MC\xa]h^\\i>C. 



(5) 



In the leading order, the FF transition operators are 
expressed as [24| 



K-i - -^A^iTiCoMo^ 

Co = V + hiiwo 

V ^ + ^M^' 
Mi ^ -gA^S < f\\ir[Ci x a]°t-\\i > C 

= -gA<f\\Tot-\\i>C, (3) 

where ^ — aZ/2R with a the fine structure constant 
and R is the nuclear charge radius, 71 = -^/l — (aZ)^, 
Cl = A/47r/(2L-h 1)11, C = l/y/2J, + 1. The prime in 
M(f ' indicates that the effects of the finite nuclear charge 



CM = ^;^l<.fl|0(A-)||z>p 

0{Q-)^gA[—+£.iS-r^t^ 
m 

0(1") = [gv— " ^{gA^ X f - igvr)]t~ 
m 

0{2-), - *^[^ X r^lVW+^J- (6) 

where gv is the vector coupling constant. Matrix ele- 
ments of the first and the second parts of O(0~), the 
first, second and the third terms of 0(1~), and 0{2~) 
correspond to Mff and Mq', £^'y, u' and x', and z in Ref. 
[T9| . respectively. The transition strengths constructed 
from the matrix elements of O(O^) and 0(1") in Eq. (6) 
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are dominant parts of Kq, while that of 0(2 ) in Eq. (6) 
becomes equivalent to Eq. (5) when A2 equals unity. 

A shell model study was made for FF /3 decays in the 
lead region, A =205~212, where an enhancement of the 
rank-zero matrix element of 75, Mq , has been found [25l |. 
The quenching of the rank-one and rank-two components 
of the decay rate due to the core polarization effects 
has been also investigated [2^ [2^. Here, the transition 
matrix elements and rates are evaluated by using Eqs. 
(l)-(5) following Refs. 

The isotones with proton number of Z =64^73, that is, 
nuclei with Uh =18^9 proton holes are considered here 
for the shell model calculations ^28^ . A closed A^ =126 
shell configuration is assumed for the parent nucleus. 
Proton holes in Ohn, and 2si/2 orbits are taken 

into account for the shell model calculations. Configura- 
tions with 2~4 and up to 2 holes are taken for 1^3/2 and 
2si/2 orbit, respectively. For the 0/iii/2 orbit, (n/i— 6) ~ 
(n/i— 4) (~12 in case of nh=\7) hole configurations are 
taken into account for <17. In case of Uh =18 {Z 
=64), 10^12 hole configurations for the O/111/2 orbit and 
additional 0^2 hole configurations for the lc?5/2 orbit are 
considered. For neutrons, the Oft.9/2, 1/5/2,7/2, 2^1/2,3/2 
and Oii3/2 orbits outside the N =82 core are taken as 
the model space, and in the /3-decays a neutron in these 
orbits changes into a proton whose orbit has holes. The 
FF transitions are induced by ;^0ii3/2 ''rO/iii/2 and 
^{fp) ~^ 7r(sd) transitions while the GT transition is 
dominantly induced by 1^0/19/2 — 7r0/iii/2 transition. In 
the /^-decays, important contributions come from transi- 
tions to low- lying states. The giant resonance region is 
energetically off the /3-decay windows even for the very 
neutron-rich cases. Though the calculations have been 
carried out with a restricted configurations, important 
states with low excitation energies are considered to be 
well described by the present shell model calculations. 

The GT transition strengths, B{GT) = Kq (eq. (2)), 
obtained by the shell model calculations arc shown in 
Fig. 1. Here, qa is taken to be quenched by a factor 
of 0.7 [H,!!^, that is, QA ^0.7xg{{'"' with g{'''"'^-1.26. 
Energies denoted are those from the parent state. As 
the proton hole number increases, the energy difference 
between the parent state and the daughter states at the 
peak of the strength and the phase space for the transi- 
tion becomes larger. Note that the transition strength is 
approximately proportional to the fifth power of the en- 
ergy difference. The summed GT strength also becomes 
generally larger for more proton holes; the sum of the cal- 
culated B{GT) values are 14.4, 14.6, 11.7, 8.5 and 5.6 for 
Z = 64, 66, 68, 70 and 72, respectively. Both of the ef- 
fects lead to shorter half-lives for the isotones with larger 
proton hole number in the /3-decays. 

We find that the most important contributions from 
the FF transitions come from the case of J'^ = 1^. We 
discuss spin-dipole strengths which are important spin- 
dependent contributions in Kq as well as the total shape 
factors averaged over electron energy. Calculated spin- 
dipoloe strengths 



B{SDX) = < X > p (7) 

with A=l are shown in Fig. 2 (upper figures). General 
feature of the strength distributions is similar to the case 
of the GT transitions. Energy difference between the 
parent state and the daughter states at the peak of the 
strength as well as the sum of the strength get larger 
for nuclei with more proton holes. The summed B{SD1) 
values are 55.5, 49.2, 48.4, 40.1 and 35.1 fm^ for Z =64, 
66, 68, 70 and 72, respectively. 

In case of FF transitions, as the shape factor C{w) de- 
pends on electron energy, it is more appropriate to discuss 
the averaged shape factor instead of Kq- Following Ref. 
(25} . the averaged shape factor is defined as 



C{w) = f/fo (8) 

with 

/Wo 
F{Z, ■w)pw{wo - wYdw, (9) 

and for FF transitions 

— - 9195 X 10^ ,^ 2x 

C{w) = — (fm^). (10) 

Averaged shape factors for 1~ transitions are shown in 
Fig. 2 (lower figures). Here the matrix elements are mul- 
tiplied by the electron Compton wave length so that they 
have dimension of 'fm'. The transitions have components 
from the hadronic vector current also, that is, the contri- 
butions from the electric dipole {El) transitions. We see 
from Fig. 2 that the the total contributions have similar 
strength distributions as those of the spin-dipole part, 
while the contributions from the spin independent part 
are large. The peak position of the transition strength 
in nuclei near Z 70 is located at /3-decay energies of 
about 6 MeV (4 MeV) for the 1" (GT) transition, which 
is found to be similar to the case of the CQRPA calcula- 
tion in Ref. [H. 

B. Half-lives of the Waiting-point Nuclei 

Decay rates and relevent matrix elements for the 
GT and FF transitions are evaluated [l^ by includ- 
ing the quenching of the axial vector coupling constant: 
.9a/3^"^'^=0.7 is taken in the present study for both the 
GT and FF transitions except for the 0~ case. In relation 
to this, similar order of quenching is found for the spin 
g factor, g^f^ =0.64 gs, in the study of spin-dipole M2 
transitions in heavy nuclei [23|. We, therefore, assume 
the same quenching in the spin-dipole transitions for 1~ 
and 2~ cases as in the GT transitions. 

As for the 0~ case, gA in {a ■ p term from 75) is 
enhanced due to the meson exchange current effects [23] . 
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FIG. 1. (Color online) GT strengths for the isotone with Z=64, 68 and 72 obtained by the shell model calculation at energies 
from the parent state. 





FIG. 2. (Color online) The same as in Fig. 1 for the spin-dipole strengths with =1 (upper figures) and the corresponding 
average shape factors (lower figures) for the isotones. 



The enhancement factor is taken to be e =gA/ 9^'^'^ =2.0 
while it is taken to be 1.0 in Mq ((j-r term) following the 
analysis in Ref. [Hi, where e =2.01±0.05 and 0.97±0.06 
are obtained for Mq and M(f , respectively, by fitting to 
the experimental /3-decay data for A =205 ^212. 

Calculated half-lives for the /3-decays of the isotones 
are shown in Table I and Fig. 3(a). Calculated /3- 
decay Q-values obtained in the shell model calculations 



are used. Calculated half-lives obtained with an approx- 
imation of using /ii =\2 =1 are also shown in Table I. 
This approximation changes the half-lives by only within 
0.3%. (The effects on the FF transition rates are within 
0.5%.) As the approximation proves to be quite accu- 
rate even for the present large Z cases, we adopt =A2 
=1 hereafter. The validity of this approximation is also 
pointed out in Ref. [25] for non-unique FF transitions. 
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TABLE I. Calculated half-lives (in units of ms) for the /J-decays of A'^ =126 isotones obtained by shell model calculations with 
the use of the quenching factor of QA/gl^'^'' =0.7 for 1+, 1~ and 2~ transitions. The enhancement factor of e =2.0 is used for 
the matrix element Mq in 0~ transitions. Results for the GT transitions as well as those for both the GT and FF transitions 
are shown for Z = 64~73. Results for the case of the approximation, ^\ =A2 =1, in the FF transitions are also shown. 
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64 


65 


66 


67 


68 


69 


70 


71 


72 


73 


GT (ms) 

GT + FF (ms) 

GT + FF (A2=l, Mi=l) (ms) 


5.76 
4.02 
4.01 


7.69 
5.31 
5.31 


11.26 
7.75 
7.74 


17.46 
10.94 
10.94 


29.31 
18.11 
18.08 


54.77 
29.49 
29.48 


102.49 
44.18 
44.07 


223.26 
84.81 
84.73 


504.79 
129.65 
129.24 


1584.4 
278.88 
278.35 




64 66 68 70 72 7^ 64 66 68 70 72 74 

z z 



FIG. 3. (Color online) (a) Calculated half-lives for the A^=126 isotones. Results of the present shell model calculations with 
GT and with GT-I-FF transitions are denoted by dashed and solid curves, respectively. The quenching factor of QA/g^'"' =0.7 
is used for both the GT and FF transitions except for 0~ transitions ( see text). Half-lives of Ref. denoted as FRDM are 
shown by a dotted curve, (b) The same as in Fig. 3(a) for the shell model calculations with contributions from both the GT 
and FF transitions. The solid curve and the curve denoted as FRDM are the same as in (a) , while the long-dashed and dotted 
curves are obtained by using a different quenching factor of Qa/q^'^'^ =0.34 with further quenching of gv / gy'^"'^ =0.67. In case 
of the dotted curve, the Q-values are increased by 1 MeV for the isotones in the FF transitions. 



Calculated half-lives obtained by the GT contributions 
only are found to be close to those in Refs. [13, [IE HH 
within ~ 10% except for Z =71 and 73. In case of Z =71 
and 73, the half-lives obtained here are shorter by about 
1.3-1.5 than those in Refs. [13, The sheU model 



interaction used here[l7[ is not the same as that in Refs. 
[lOl [l6j . The half- lives become shorter for isotones with 
smaller Z as both the strength and the energy difference 
between initial and final states get larger. The FF tran- 
sitions are found to be important to reduce the half-lives 
by 1.4—5.7 times of those given by the GT contributions 
only. Their contributions become more important com- 
pared to the GT ones for larger Z cases. The dominant 
contributions come from the 1" transitions. The con- 
tributions from the 0^ transitions are about 10% of the 
total FF contributions; 7-8 % for Z =64-68, 9-10 % for 
Z =69-71 and 11% (13%) for Z =72 (73). The reduction 
of the half- lives due to the inclusion of the FF transitions 
is not as large as that obtained in Ref. [l^. The half- 
lives obtained here are, however, short compared with 
the standard data of Ref. ^12] except for Z =71 usu- 



ally employed in nucleosynthesis network calculations as 
shown in Fig. 3(a). The present half-lives of the shell 
model calculations are shorter than those of the stan- 
dard values by 2.3—8.3 for even Z and by 1.4—2.0 for 
odd Z (except for Z =71), respectively. They increase 
monotonically as Z increases showing no odd-even stag- 
gering found in FRDM's. This is due to the absence of 
the odd-even staggering in the /?-decay Q-values in the 
shell model calculations in contrast to the case of FRDM 
calculations. 



The present half- lives are longer than those of Ref. [15 1 
by about 1.1-1.3 (1.5) for Z =64-67 (68) and by twice 
for Z =69 and 70, respectively. They are, on the other 
hand, short compared to the half- lives of Ref. [Slj , where 
similar CQRPA calculations have been done as in Ref. 
[l5| but without energy dependent smearing of the GT 
and FF transition strengths. The present half-lives are 
shorter than those of Ref. 3l[ by about a factor of 1.5 
for Z =69, 70 and about by twice for Z = 72, 73, while 
they are close to each other for Z =64—68 and 71. As the 
shell-model calculations include the spreading of the GT 
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TABLE II. Values of log fot and y C{w) for FF transitions in nuclei with Z =79, 80 and A'^ =126. Calculated values obtained 
by using quenching factors of {gA/gA^'^, gv/gy''^'^) for 1~ transitions and enhancement factor e for 0^ transitions as well as 

experimental valuesjs^] are shown. Values of Ref. ^] are also shown for Z =80. Values of log fot and \J C{w) in the parentheses 
are for the case in which with proton holes are restricted to Ohii/2, ld5/2 and 2si/2 orbits. 



Transitions 


Values of gA and gv 


log fot 


^C(w) (fm) 


Initial Final 


e for 0" 


(gA/gi 


ree 


/ tree\ 

gv/gv ) 


for 1" 


^"•"Hg (0+) ^""Tl (0- g.s.) 

Expt. [32] 


2.0 
1.8 
Ref. [25] 




5.199 (5.087; 
5.432 (5.320; 
5.173 
5.41 




76.3 (86.8; 
58.3 (66.3; 

78.6 

59.8 




^"''Hg (0+) ^""Tl (1", 0.3049 Me V) 
Expt. [32] 






% 


,0.34, 0.67; 
(0.51, 0.30 
0.47, 0.64) 
34, 0.40) 
Ref. [25] 




5.017 (4.929; 
5.157 (5.127; 
4.921 (4.832; 
5.267 (5.178; 

5.181 

5.24 




94.0 (104.1) 
80.0 (82.9) 
105.0 (116.4) 
70.5 (78.2) 
77.9 
72.7 


^"=Au (3/2+) ""^'Hg (1/2-) 
Expt. [32.1 






[a) (0.34, 0.67; 
b) (0.51, 0.30 
(c) (0.47, 0.64) 
(0.34, 0.94) 




6.197 
8.171 
6.412 
5.793 
5. 


6.116; 
7.834; 
6.326; 
5.726; 
79 




24.2 
2.49 
18.9 
38.7 

38 


26.5; 
3.67; 
20.8; 
41.6; 
.6 




''"''Au (3/2+) ^"''Hg (3/2j", 0.4675 MeV) 
0" 

o~ + 1- 

Expt.[32l 


2.0 
1.45 
1.45 
1.45 
1.45 
1.45 




'a) (0.34, 0.67; 
b) (0.51, 0.30 
(c) (0.47, 0.64) 
(0.34, 0.94) 




5.674 (5.54i; 
6.832 (6.699; 
6.502 (6.247; 
6.000 (5.870; 
6.266 (6.073; 
6.425 (6.173; 
6.43 




44.1 (51.5^ 
11.6 (13.6; 
17.0 (22.8^ 
30.3 (35.2; 
22.3 (27.9; 
18.6 (24.9; 
18.5 




^"^Au (3/2+) ^"''Hg (5/2j-, 0.3792 MeV) 
Expt. [32] 




(a) (0.34, 0.67) 

(b) (0.51, 0.30) 

(c) (0.47, 0.64) 
(0.10, 0.15) 


5.108 (5.006; 
5.409 (5.308; 
5.057 (4.956 
6.343 (6.242; 
6.37 




84.7 (95.2) 
59.9 (67.3) 

89.8 (100.9) 
20.4 (23.0) 

19.8 


^"''Au (3/2+) ^"''Hg (3/2^7, 1.2806 MeV) 

0- + 1- 

Expt. [32] 


1.5 
1.5 
1.5 
1.5 


(a) (0.34, 0.67) 

(b) (0.51, 0.30) 

(c) (0.47, 0.64) 
(0.34, 0.94) 


5.497 (5.457; 
5.556 (5.454; 
5.433 (5.413 
5.364 (5.374; 
5.51 




54.1 (56.7) 
50.6 (56.8; 
58.3 (59.6 
63.0 (62.4; 
53.3 




^'■"'Au (3/2+) ^'■"'Hg (5/2^, 1.2806 MeV) 
Expt. [32] 








'0.34, 0.67; 
(0.51, 0.30 
0.47, 0.64) 
34, 0.94) 




5.397 
5.646 
5.333 
5.176 
5. 


5.8O7; 
6.064; 
5.746; 
5.584; 
51 




60.7 (37.9; 
45.6 (28.2; 
65.4 (40.6; 
78.3 (48.9; 
53.3 




^'^^Au (3/2+) ^"''Hg (1/2^, 1.4472 MeV) 
Expt. [321 






'1o 


(0.34, 0.67; 
(0.51, 0.30 
0.47, 0.64) 
34, 0.94) 




6.632 (6.163; 
7.191 (7.956; 
6.708 (6.386; 
6.300 (5.767; 
6.29 




14.7 (25.1; 
7.7 (3.2) 

13.4 (19.4; 

21.5 (39.7; 
21.7 





and FF transition strengths, it is reasonable that shorter 
half-lives are obtained in the present calculations. 

We here consider other possible quenching factors for 
1~ transitions. A large quenching of qa and gv for 
1~ spin-dipole transitions was suggested by studies of 
FF /3-d ecay s and related processes in the lead region 
[11 [ig. In Ref. [H, effective quenching factors 

of Qa/Qa'^'^ ^0.47 and gv / Qv'^^ ^0.64 due to core po- 
larization effects are obtained for 1^ transitions. In 
Ref. [18], two sets of quenching factors, qa/ Qa'^'^ =0.34, 
gvlgv'^'^O.Ql and qaIq^a'"' =0.51, gv/gv"" =0.30, are 
obtained from the analysis of the FF transition, ^osq^j 



(1/2+, g.s.) ^ 205pb (1/2-). 

We study FF /3-decays in nuclei with Z =78~80 and 
N =126, where ex per imental data are available for 0~ 
and 1~ transitions [321 ■ 



Calculated values of log f^t and yC{w) (fm) for the 
FF ^-decays, ^^^Hg ^o^Tl (0", 1~) and ^o^Au ^osjjg 
(1/2", 3/2", 5/2"), are given in Table II for the follow- 
ing sets of the quenching factors; (qa I Q^I^'^ , gv/gv'^'^) = 
(a) (0.34, 0.67) fl^, (b) (0.51, 0.30) [3 and (c) (0.47, 
0.64) [2^ for 1~ transitions. Experimental values jl^l 
as well as the calculated values given in Ref. [25} for Z 
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TABLE III. Half-life for the /3-decay of ^"^Pt. Calculated values obtained by using quenching factors of {gA/g^^'^ , gv /gy'^") for 
1~ transitions and enhancement factor e for 0~ transitions as well as experimental value [33] are shown. Values of the half-life 
in the parentheses are for the case in which proton holes are restricted to Ohii/2, ld5/2 and 2si/2 orbits. . 



Transitions 


Values of gA and gv 


Half-life (s) 


Initial Final 


e for 0" 


{9A/gA^% gv/gv"^") for 1 


"""Pt (0+) ""*Au (0-) 


2.0 




65.1 (62.0) 


^""Pt (0+) ^"''Au (0- + 1-) 


2.0 


(a) (0.34, 0.67) 


22.8 (18.6) 




2.0 


(b) (0.51, 0.30) 


31.8 (26.9) 




2.0 


(c) (0.47, 0.64) 


21.1 (17.1) 




2.0 


(0.7, 1.0) 


10.9 (8.6) 




2.0 


(0.34, 1.0) 


14.4 (11.3) 




2.0 


(0.51, 1.0) 


12.7 (10.0) 


Expt.[32l 






10.3±1.4 



=80 and cases for some other sets of quenching such as 
(0.34, 0.40) and (0.34, 0.94) are also given. In case of 
0~ transitions, the enhancement factor of e = gA/gA^'^'^ 
=1.45^2.0 is adopted for the matrix element, . Val- 
ues in the parentheses in the Table are for the case in 
which proton holes are restricted to Ohxj^/2: 1^3/2 and 
2si/2 orbits. 

In case of Z =80, for the 0^ transition, calculated log 

/oi and \J C{w) values for e =2.0 are close to those of 
Ref. and become closer to the experimental values 
for e =1.8. For the 1~ transition, the calculated log f^t 
values for the sets (a), (b) and (c) are smaller than those 
of Ref. [13 by 0.17, 0.02 and 0.26, respectively. They get 
closer to the observed values if a set with more quenching 
for gv, (0.34, 0.40), is used. 

In case of Z=79, transitions to the l/2j^ and 
states in ^^^Hg are dominantly induced by 1^ transitions, 
while the transition to the 3/2^ state is a mixture of 0^ 
and 1~ transitions. For the 1~ transition to the 1/2^ 
(g.s.) state, the three sets of the quenching give larger 

(smaller) log f^t {\J C{w)) values than the experimental 
ones. A set with smaller quenching for gy, (0.34, 0.94), 
gives the values close to the experiment. For the transi- 
tion to the 3/2;," state (E^ =0.4675 MeV), the 0" transi- 
tion proceeds too fast for e =2.0. When an enhancement 
factor of e =1.45 is used, the three sets of the quenching 
for the 1^ transitions can give log f^t and ^J-^^ values 

close to the experimental ones. The set used for the 1/2 
case, (0.34, 0.94), works also in this transition. As for the 
5/2 J" state, any of the three sets of the quenching for the 
1~ transition give transitions too fast compared to the 
experiment. In this case, an exceptionally large quench- 
ing, (0.10, 0.15), is necessary to explain the experiment. 

Although the spins of the excited states of ^'^^Hg at 
=1.2806 MeV and E^ =1.4472 MeV are not deter- 
mined yet, the present analysis suggests that the state 
at Ex =1.2806 MeV is 3/2" or 5/2" and the state at 
E.^ =1.4472 MeV is 1/2" (see Table II). Other possibil- 
ities, that is, 1/2^ (3/2^ and 5/2^) for the state at E^ 
=1.2806 MeV (1.4472 MeV) can be excluded. 

In case of Z =78, the half-life for the transition, ^""'Pt 



(0+) —7' ^'''^Au has been obtained by measuring two 7 
decays in Au The averaged experimental half-life 
is 10.3 s. Here, the half-life is calculated by taking into 
account the FF 0~ and 1~ transitions. The enhancement 
factor of e =2.0 is used for the 0~ transition. The half- 
lives obtained with the three sets of the quenching for 
the 1~ transitions are given in Table III. The calculated 
half-lives are longer than the observed one by about 2^3. 
The values become close to the observed one for the cases 
without quenching for gy (see Table III). The present 
analysis suggests smaller quenching for gv for Z =78. 

The present analysis for Z =78~80 suggests a large 
quenching for gA- A large quenching for gy is not nec- 
essarily needed for Z =78 and for most cases oi Z =79. 
We study dependence of the calculated half-lives on the 
quenching of gA and gv ■ Half- lives obtained with the set 
(a) for the 1~ transitions are shown in Fig. 3(b). The en- 
hancement factor e for the 0^ transitions is not changed; 
e =2.0. The half-lives get longer but still remain shorter 
than those of Ref. [l^l except for Z =71. They become 
close to the values of Ref. [3l[ . 

Although the experimental Q-values are used for the 
analysis of the /3-decays for Z =78^80, the present shell- 
model calculations give smaller Q- values compared to the 
experimental ones. The differences are 1.3, 1.6 and 0.9 
MeV for Z =80, 79 and 78, respectively. As there are no 
experimental information on the Q- values for Z =64^-^73, 
we used the calculated Q- values. They are small com- 
pared to the mass differences between the parent and 
daughter nuclei in Ref. [Tl| by about 2 MeV (3 MeV) 
for even (odd) Z nuclei. Now, we study the dependence 
of the half-lives on the Q-values of the transitions. Cal- 
culated half-lives are shown in Fig. 3(b) for the case, 
in which the Q-values obtained by the shell model cal- 
culations are increased by 1 MeV for the isotones with 
Z =64^^73. The effects of the change of the Q-values, 
which results in the enhancement of the phase space vol- 
umes, are found to be large. The half-lives get shorter 
than the case in Fig. 3(a) with the quenching, (0.7, 1.0) , 
but they do not become as short as those of Ref. [iJl ■ 
They are still longer than those of Ref. [l^ by about 
30% for Z =66-68 and 50% (60%) for Z =69 (70). It is 
quite important to obtain experimental information on 
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the masses of waiting point nuclei and their daughters. 



III. R-PROCESS NUCLEOSYNTHESIS 

We discuss possible effects of the short half-lives of the 
waiting point nuclei at 7V=126 obtained in the preced- 
ing section on the r-process nucleosynthesis. The depen- 
dence of the abundances of the elements around mass 
number A~195 on the half-lives of the nuclei is inves- 
tigated for various astrophysical conditions. We use an 
analytic model for neutrino-driven winds [ssl ] for the time 
evolution of thermal profiles. The wind solutions which 
pass the sonic points are adopted. We take into account 
the relation of the neutrino luminosity L^j and the mean 
neutrino energy to the entropy per baryon S/k and the 
mass ejection rate M of the winds in accordance with 
[33 |. The neutrino energy spectra are assumed to obey 
Fermi distributions with zero chemical potentials. The 
temperatures of Vg, 9^, and v^: = v^^t and v^^r are set 
to be (r^^,Tp^,r^J = (3.2 MeV, 5 MeV, 6 MeV) [H. 
The luminosity of each fiavor of neutrinos is equiparti- 
tioned and is taken to be — (0.5 ^ 5.0) x 10^^ erg 



s . The neutrino luminosity, entropy per baryon, and 
mass ejection rate of each wind model are listed in Table 
IV. 

We calculate r-process nucleosynthesis using a reaction 
network consisting of 3517 species of nuclei. The initial 
electron to baryon ratio is taken to be Yg =0.40. Neutrino 
processes on nucleon and "^He are included in the net- 
work calculations. Charged current reactions, n(i'e, e~)p, 
p(i/e, e"'")n, ''He (i/e,e~p) ^He, and '^He (j/gje+n) '^H give 
dominant contributions. Besides them, charged current 
reactions, ''He (I'e, e~pp) ^H, ''He (i/e, e+nn) ^H, and neu- 
tral current reactions, ''He [v, v'x\) ■^He, ''He (i^, i^'p) '^H, 
''He (i^, v'A) ^H, ''He {v, z/'nnp) ^H and those induced by 
i/'s are included. Cross sections of ref. [s^ obtained by 
the WBP Hamiltonian [s^] are used. The cross sections 
of n and p are adopted from [s^ . The neutrino processes 
suppress the production of the elements around the third 
peak in the r-process as neutrons are changed into pro- 
tons, which results in the increase of Ke, and the produc- 
tion of seed elements is enhanced through the neutrino- 
induced breakup of ''He [9]. We artificially shorten the 
time scale of the explosion by a constant factor /( in order 
to proceed the r-process to the synthesis of the 3rd peak 
elements because it was found that the shorter expansion 
time scale tends to decrease the roles of neutrino inter- 
actions and keep low Ye in favor of the r-process We 
set the final temperature of the winds to be = 8 x 10* 
K. The duration time r where the temperature decreases 
from 5 X 10^ K to 2 X 10^ K is T = (1.8 - 5.6) ms. 
The multiplying factor ft and the duration time t for 
r-process calculation are also shown in Table IV. 

The abundance of the elements obtained by using the 
present short half-lives in the r-process element network 
are compared with those obtained with the use of the 
standard data of ref. [T3|. We show in Fig. 4 a case 



in which the ratio of the height of the third peak over 
that of the second peak in the r-process is close to the 
solar abundance ratio, 3~4. Parameters chosen are Lj/,5i 
=0.5 where L^^si =LJ{IQ^^ erg s'^), S/k = 133 and 
T = 5.6 ms. Here, the half-lives obtained for the case 



with the quenching of Qa/Qa^'^ =0.7 (gv/gy'^'^ =1) for 
1~ transitions are used. Other inputs except for the half- 
lives are not changed. We show in Fig. 5 the calculated 
results for the case of a different set of the quenching 
factors for 1~ transitions, (g^i/ff^'^^'^, 5v/.9v^^^) =(0.34, 
0.67). A case with Q-values increased by 1 MeV is also 
shown. We find a slight shift of the third peak of the ele- 
ment abundances toward higher mass region. This shows 
more rapid build-up of the r-process elements caused by 
the shorter half-lives of the waiting point nuclei. The iso- 
tones with higher Z accumulate more abundantly at the 
freezing time of the neutron-capture flow, which by suc- 
cessive /3-decays naturally results in larger abundances of 
the elements at higher mass region. 

When the decrease of the temperature becomes slower 
or the entropy of the system is larger, the magnitude of 
the shift of the peak gets larger. The shift of the mass 
number at the third peak region is defined by 



AA^ <A >Mod - < A >STD 

j:aA-y{a) 



<A> 



(11) 



where Y{A) is the abundance of the element with 
mass number A and the summation is taken over 189< 
A <203. 'Mod' and 'STD' refer to the cases of modi- 
fied half- lives and standard ones, respectively. The mass 
shifts at the third peak of the r-process are shown in Ta- 
ble IV for several diff'erent values of L^^^i, S/k and r. 
Here, the condition is being kept so that the ratio of the 
height of the third peak over that of the second peak is 
close to the solar abundances ratio. We see from Table 
IV that the mass sfift AA is about unity in all cases. 
Dependence on the astrophysical condition is found to 
be small. Dependence on the difference of the quenching 
factors is also as small as 25~30 %. The decreace of the 
mass sfift due to larger quenchings is recovered or even 
reversed to an increase by the change of the Q-values. We 
also find that the mass shift at the second peak, where 
the average is taken over 124< A <136, is quite small; 
AA = -0.006 0.002. We thus emphasize that the shift 
of the third peak of element abundance toward higher 
mass region by about unity is a common phenomenon in 
the present astrophysical conditions. When half-lives of 
the N =82 isotones are changed to the values obtained 
by shell model calculations [ll[ , effects on the shift of the 
third peak is found to be quite small. The mass shift AA 
decreases only by 6% while the mass shift at the second 
peak is AA =0.26^0.28. The magnitude of the shift at 
the third peak is rather modest, but it is a robust effect 
and is of the same order as that caused by variations 
of mass formulae. As we see from Table IV, the posi- 
tion of the third peak of the solar abundance, A ^190, 
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TABLE IV. The mass shift at the third peak of the element abundance in the r-process defined by eq. (8) as well as the peak 
position for the case of the modified half-lives are shown for several astrophysical conditions. The ratio (R) of the height of the 
second peak over that of the third peak is also given. For each astrophysical condition, the first row corresponds to the case 
of the quenching of (<7a/<7a'^'^'^, gv/gy'"'"^) =(0.70, 0.0), while the second and the third rows correspond to the case of (a) (0.34, 
0.67) for 1~ tarnsitions. The Q-value is increased by 1 MeV for the third row (denoted as (a)*). 





S/k 






r (ms) 


ratio R < A >Mod A A 


0.5 


133.4 


2.34 X 10"" 


0.06 


5.60 


3.45 196.972 1.101 
(a) 3.21 196.684 0.813 
(a)* 3.84 197.218 1.347 


1.0 


118.9 


7.43 X 10"" 


0.08 


4.04 


4.23 197.118 1.058 
(a) 3.89 196.837 0.777 
(a)* 4.75 197.359 1.299 


2.0 


105.9 


2.36 X lO"'' 


0.10 


2.90 


4.77 197.234 1.008 
(a) 4.35 196.965 0.739 
(a)* 5.03 197.466 1.240 


5.0 


90.91 


1.08 X 10"* 


0.11 


1.78 


3.80 197.313 0.949 
(a) 3.52 197.066 0.702 
(a)* 3.94 197.523 1.159 




100 120 140 160 180 200 180 185 190 195 200 205 210 

Mass Number A Mass Number A 

FIG. 4. The abundances of elements in the r-process nucleosynthesis obtained by using the present /3-decay half-lives for the 
A''=126 isotones (denoted as 'modified') and standard half-lives of ref. [l2| (denoted as 'standard'). The quenching factor of 
gA/gy^'^" = 0.7 {gv / gv'^'' =1.0) is used both for the GT and FF (1" and 2~) transitions. The right figure displays the same 
result in the restricted range of A >180. 



is not necessarily reproduced by the present calculations. 
However, we should note that the third peak region can 
be affected by various effects besides the half-lives of the 
/3-decays as we pointed out. Neutron-star merger and 
gamma-ray bursts are considered to be possible alterna- 
tive r-process sites [1, |^ . The assignment of the site of r- 
process including those besides the supernova explosions 
is an important and interesting subject to be explored 
further in future. 



IV. SUMMARY 

In summary, we have evaluated half-lives of the 
/3-decays of N =126 isotones taking into account both 
the GT and FF transitions. The half-lives have been 
found to be reduced by including the FF transitions. 
The effects of the short half-lives obtained here on the 



nucleosynthsis of the r-process are investigated. The 
third peak of the abundance of the elements in the r- 
process has been found to be shifted toward higher mass 
region. Although the magnitude of the shift is rather 
modest, it is found to be a robust effect independent of 
the present astrophysical conditions for the r-process as 
well as the quenching factors of gA and gv adopted in 
the shell model calculations. The results obtained here 
might be improved quantitatively by overcoming the 
limitations in the present structure calculations such as 
configuration space. Improvements of various nuclear 
properties for the input for the nucleosynthesis network 
of the r-process in addition to the /3-decay half-lives, 
especially the masses of nuclei at and around the waiting 
points, are important issues to be made in future studies. 

The authors would like to thank Y. Utsuno for use- 
ful discussions. This work has been supported in part 
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180 185 190 195 200 205 210 180 185 190 195 200 205 210 

Mass Number A Mass Number A 

FIG. 5. The same as in Fig. 4 ai A >180 for the quenching factors of {gA/g^"'', gv /gv"^^) =(0.34, 0.67) for the 1^ transitions. 
The figure on the right side includes the effect of the increase of the Q-values. 
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